Non-Gaussian entanglement distillation for continuous variables 
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Entanglement distillation is an essential ingredient for long distance quantum communications 
In the continuous variable setting, Gaussian states play major roles in quantum teleportation, 
quantum cloning and quantum cryptography !2|. However, entanglement distillation from Gaussian 
states has not yet been demonstrated. It is made difficult by the no-go theorem stating that no 
Gaussian operation can distill Gaussian states [3, 0, 0]. Here we demonstrate the entanglement 
distillation from Gaussian states by using measurement-induced non-Gaussian operations, circum- 
venting the fundamental restriction of the no-go theorem. We observed a gain of entanglement as a 
result of conditional local subtraction of a single photon or two photons from a two-mode Gaussian 
state. Furthermore we confirmed that two-photon subtraction also improves Gaussian-like entangle- 
ment as specified by the Einstein-Podolsky-Rosen (EPR) correlation. This distilled entanglement 
can be further employed to downstream applications such as high fidelity quantum teleportation [f| 
and a loophole- free Bell test 0]. 
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Long distance quantum communications rely on the 
ability to faithfully distribute entanglement between dis- 
tant locations. However, inevitable decoherence and the 
inability to amplify quantum signals hinder efforts to ex- 
tend a quantum optical link to the practically large scale. 
To overcome this difficulty, entanglement distillation can 
be used - a protocol in which each distant party locally 
manipulates particles of less entangled pairs with the aid 
of classical communication to extract a smaller number 
of pairs of higher entanglement [8| . In discrete variable 
systems many distillation experiments have already been 
demonstrated [IfMCklEl- 

An alternative to the discrete variable system is the one 
described by continuous variables (CV), typically repre- 
sented by the quadrature amplitudes of a light field. For 
CV quantum information, Gaussian states and Gaussian 
operations [13jj are of particular importance. They are 
readily available in the laboratory and serve as a com- 
plete framework for many quantum protocols Q. The 
first experiments of CV entanglement distillation were 
reported in [3, EH • These two works rely on Gaussian 
operations. However, it was theoretically proven that 
Gaussian operations can never distill entanglement from 
Gaussian state inputs - this is known as the no-go the- 
orem of Gaussian operations [1, 0, Q. In [3, [HI], the 
inputs had been subjected to some specific classes of non- 
Gaussian noise, such as p hase-diffusion [Tij or temporally 
varying attenuation [15]. In those cases, well established 
Gaussian technologies can be applied to distill the entan- 
glement. 

So far, there has been no demonstration of entangle- 
ment distillation with Gaussian inputs. This task essen- 
tially requires a new technology of non-Gaussian opera- 
tions. Recent theories also revealed that this is a must 
to realize quantum speed-up of CV quantum information 
processing (QIP) [161 ] . Triggered by this new paradigm 



of non-Gaussian QIP, the research field extending to the 
non-Gaussian reg ime has rapidly developed [l7l Il8l . [l^ ] . 
The increase [20( and preparation 21( of entanglement 
from Gaussian inputs by non-local photon subtraction 
were also demonstrated. These are important steps to- 
wards the realization of entanglement distillation from 
Gaussian states. 

Here we report on the entanglement distillation di- 
rectly from CV Gaussian states by using local pho- 
ton subtraction as non-Gaussian operations, circumvent- 
ing the no-go restriction on Gaussian operations. A 
schematic of our experiment is depicted in Fig. [T] A 
continuous wave squeezed vacuum is generated from an 
optical parametric oscillator (OPO) detailed elsewhere 
[191 ] . The initial Gaussian entangled state is prepared by 
splitting the squeezed vacuum by half at the first beam 
splitter and is distributed to the separate parties, Alice 
and Bob. This half-split squeezed vacuum with squeez- 
ing parameter r is effectively equivalent to the two-mode 
squeezed vacuum with r/2 (they are compatible by lo- 
cal unitary operations, see Appendix I A 2[) . At each site 
of Alice and Bob, a probabilistic non-Gaussian opera- 
tion - photon subtraction - is performed. Specifically, 
a small part of the beam is picked off by a polarizing 
beam splitter with the variable reflectance R and sent 
through filtering cavities [l9[ to an avalanche photodi- 
ode (APD) to detect a photon (Fig. [TJ. Each photon 
detection at the APD heralds a local success of the pho- 
ton subtraction attempt. Conditioned on the subtraction 
of a photon by a single party (single-photon subtraction) 
or the simultaneous subtraction of a photon by both par- 
ties (two-photon subtraction) , Alice and Bob retain those 
two-mode states which have successfully had their entan- 
glement increased. While the single-photon subtraction 
scheme will have a higher success rate, the two-photon 
subtraction scheme will give a more Gaussian-like final 
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a decomposed way as 

W(x a ,Pa,x b ,Pb) = W s {x-,p-)W v (x + ,p + ), (1) 
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FIG. 1: Schematic of the experiment. OPO: optical para- 
metric oscillator, BS: beam splitter, HBS: half-beam split- 
ter, BHD: balanced homodyne detector, LO: local oscillator, 
FC: filtering cavity, APD: avalanche photo diode, PZT: Piezo 
electric transducer. The phases of the LOs, 9a and 8b are 
stabilized by electronic feedback to the PZTs through field 
programmable gate array (FPGA) modules. The outputs of 
the BHDs are recorded by a digital oscilloscope triggered by 
either logical "AND" or "OR" of the click signals from the 
two APDs. For the state verification, a set of the homodyne 
outcomes are numerically converted into the "+/— " basis. 



state which is more readily applicable to further process- 
ing such as e.g. quantum teleportation. 

The distillation works since the local photon subtrac- 
tion changes the non-local unfactorizable correlations of 
the initial state. To see this intuitively, let us describe the 
initial squeezed state with squeezing parameter r in pho- 
ton number basis as |0,0)-^-(|0,2)+-\/2|l, 1) + |2,0)) + 

0(A 2 ) where \m, n) — \m) A\n)B, A = tanhr and we omit- 
ted the normalization. For small A, it is almost factoriz- 
able since |0, 0) is dominant. Applying the single-photon 
subtraction, represented by an annihilation operator a, 
the state is transformed to be \/2(|0, 1) + |1, 0)) + 0(A) 
which is clearly more entangled - the first Bell state term 
corresponds to 1 ebit of entanglement. We emphasize 
that the term O(A) is not negligible for larger A and crit- 
ically contributes to the distillation, and thus the scheme 
works for any A in principle. See Appendix IA 31 for the 
rigorous formulation and for the two-photon subtraction. 

The verification of the distilled (or undistilled) states 
is carried out by a quantum tomographic method with 
two local homodyne measurements (see Fig. [I}. For to- 
mography, in general the local oscillator (LO) phases 9 a 
and 9b (for Alice and Bob's detectors, respectively) have 
to be swept over all possible combinations to collect full 
information of the two-mode state. In our case, how- 
ever, it is significantly simplified due to the fact that our 
states, distilled or undistilled, are always represented in 



where W is the Wigner function for the two-mode state, 
xa,xb,Pa,Pb are the quadrature amplitudes of modes 
A and B, x± = £A ^ La , p± = EA ^§ M - and W s , W v are the 
Wigner functions for a (zero-, one-, or two-) photon sub- 
tracted squeezed vacuum and the vacuum respectively. 
For such a state, the scans of the homodyne measure- 
ments are necessary only for 9a = 9b and the exper- 
imental data of x± is numerically obtainable from the 
measured xa and x B (Fig. [T]) . It should be stressed that 
although we assume the state factorization of ((T|), it can 
be directly assessed by verifying experimentally whether 
the state of the g+h mode is indeed a pure vacuum state. 
For details, see Appendix I A 11 

Examples of reconstructed Wigner functions obtained 
by the single- and two-photon subtractions, as well 
as the initial squeezed state are shown in Fig. [2k- 
c. The outputs of the homodyne detectors were sam- 
pled at 6 different phases of LO, namely 9a — 9b — 
0, 7r/6, 7r/3, 7t/2, 27t/3, 5ir/6. We extract the measured 
values of the quadratures xa and xb by applying a mode 
function to the recorded traces [H, EH Elj]. After cal- 
culating the corresponding values of x± , we reconstruct 
the density matrices for the "+" and "— " modes by the 
conventional maximum likelihood estimation (23j with- 
out any correction of detection losses. 

As shown in Fig. [2fe, for the "+" mode states we got 
almost perfectly pure vacuum states with more than 99% 
accuracy. We confirmed that this holds irrespective of the 
initial squeezing level. This experimental evidence justi- 
fies our tomography scheme based on the relation fTJ). 
On the other hand, for the "— " mode states we observed 
two different kinds of non-Gaussian state depending on 
whether single photon or two photons were subtracted 
(Fig. [2k and b). They respectively correspond to the 
odd and even Schrodinger cat state, i.e. \a) — |— a) and 
| a) + |— a) where \a) is a coherent state with coherent 
amplitude a. Having these reconstructed states we can 
use the relation |T]) backwards to calculate the amount 
of entanglement shared by Alice and Bob. Specifically, 
we calculate the logarithmic negativity E^f which is a 
monotone measure of entanglement [241 ]. 

Fig. [2JI shows the experimental negativities of the 
undistilled Gaussian states, the states distilled by single- 
photon subtraction with R = 5%, and by two-photon 
subtraction with R — 10% as functions of the squeez- 
ing of the initial input states. When evaluating negativ- 
ity, one must take care of its strong dependency on the 
size of the data set. We investigated the behavior of the 
negativity on the data size and deduced an extrapola- 
tive value corresponding to an infinitely large data set 
for each point in Fig. [2JI (See Appendix iBj) . Note that 
without this analysis, evaluation of negativity with finite 
sized data very likely goes into an overestimate of the 
negativity. As shown in the figure, over a wide range of 
the initial squeezing we got clear gains of entanglement 
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FIG. 2: a, b, c Experimentally reconstructed Wigner functions and their contour plots of the "— " mode states for (a) 
Distilled state via single-photon subtraction with R =5%, (b) Distilled state via two-photon subtraction with R =10%, and 
(c) Undistilled initial state (squeezed vacuum with R =0%), all with initial squeezing of -3.2 dB. d Experimental logarithmic 
negativities as functions of the initial input squeezing. Here the initial squeezing refers to the squeezing of the states right 
after the OPO and is deduced from separately measured classical parametric amplification of the OPO. For the single-photon 
subtracted states (red) and the undistilled states (green), 600,000 samples over 6 phases were used for the reconstruction of 
each point. For the two-photon subtracted states (blue) 18,000 - 48,000 samples were used. The dashed curves are theoretical 
predictions based on independently measured experimental parameters. Every error bar represents an uncertainty of the state 
reconstruction and was estimated via a MonteCarlo simulation using the corresponding experimental parameters, e, The 
photon number distributions of the experimentally reconstructed "+" mode states corresponding to (a) and (b). 



relative to the undistilled Gaussian states for both the 
single- and two-photon subtracted schemes. 

A practical difference between the single- and two- 
photon subtracted scheme is on their rates of event detec- 
tion. In the single-photon experiment the rate is around 
a few thousands per second, but in the two-photon ex- 
periment there are only a few events per second. So while 
for the former we can use hundreds of thousands of sam- 
ples for the state reconstruction, for the latter we can 
only use a few tens of thousands limited by the long-term 
stability of the setup. In Fig. the experimental neg- 
ativities for the single-photon subtracted states and the 
undistilled states are in very good agreement with the- 
ory, but ones for the two-photon subtraction are slightly 
below the theoretical predictions. This may be due to an 
uncontrollable drift of the system during a long period of 



the measurements. 

As can be seen in Fig. [5Ji, in terms of the logarith- 
mic negativity the two-photon subtracted scheme docs 
not have an advantage over the single-photon subtracted 
scheme despite its significantly lower success rate. How- 
ever the two-photon subtracted distillation transforms 
a two-mode Gaussian state into one relatively close to 
a Gaussian state (see Fig. (5b)- Hence one would ex- 
pect that states distilled by this scheme still possess 
a Gaussian-like property of entanglement. For Gaus- 
sian states, two-mode entanglement is usually specified 
in terms of the Einstein-Podolsky-Rosen (EPR) corre- 
lation quantified by ((A£_) 2 )((Ap + ) 2 ). Since the "+" 
mode is always a vacuum state (see Eq. [TJ, we can fo- 
cus on the degree of squeezing of the "— " mode as an 
equivalent measure. We carried out measurements of the 
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FIG. 3: Squeezed variances of x- (normalized by the vacuum 
level). For the undistilled states each point was calculated 
from the full-reconstructed density matrices common to the 
negativity measurements. For the distilled states each point 
was obtained by directly measuring variance of 1,600 - 5,000 
samples at the most squeezed phase. 



variance of x— at its most squeezed phase conditioned on 
two-photon subtraction ( Fig. [3]). Note that this mea- 
surement is considerably faster than reconstruction of a 
full state and is possibly more accurate due to its sim- 
plicity. The results in Fig. [3] show that the two-photon 
subtraction improves the EPR correlation of a state with 
up to 4 dB of initial squeezing. An improvement on the 
EPR correlation gives us an operational measure of suc- 
cess of the two-photon subtracted distillation in terms 
of the fidelity of CV quantum teleportation Q . On the 
other hand, the single-photon subtracted distillation is 
never able to improve the degree of the EPR correlation. 
This fact demonstrates a clearly different nature of the 
distilled entanglement between the two schemes. 

In conclusion, we have for the first time demonstrated 
CV entanglement distillation from Gaussian states by 
conditional local photon subtraction. Because of the im- 
portance of Gaussian states and the no-go theorem of the 
Gaussian distillation, this has been a long-standing ex- 
perimental milestone to be achieved. Our scheme would 
serve as the de-Gaussifying process of a more generic 
distillation protocol proposed in (25j by combining it 
with the already demonstrated Gaussification processes 
0, EH ) which would realize long-distance CV quantum 
communications. Finally, our non-Gaussian entangled 
states are not only useful for communications but also 
for fundamental problems such as a loophole-free test of 
Bell's inequality [7j. 

H.T. acknowledges the financial support from G-COE 
program commissioned by the MEXT of Japan. 



APPENDIX A: THEORETICAL DESCRIPTION 
OF THE DISTILLED STATES 

1. Factorizability of the distilled states 

Here we discuss the local photon subtraction from a 
two-mode entangled Gaussian state made by splitting a 
squeezed vacuum. Fig. 2^, illustrates this scheme, where 
one of the inputs for the half-beam splitter is in the vac- 
uum state. Conditioning on either photon detection at 
one of the detectors or coincidence detection of two de- 
tectors heralds a non-Gaussian two-mode state effectively 
described as follows (apart from the normalization). 

Hout) AB = a n /a B B B AB ^/4)S A (r) \0) A \0) B (Al) 
= B AB (n/4) a n A A+nB S A (r) \0) A \0) B , (A2) 

where S A (r) is the squeezing operator with the squeez- 
ing parameter r and B ab (tt/4) is the half-beam split- 
ter operator for the modes A and B. More generally 
we denote the beam splitter operator with reflectance 
R as B AB {9) where 8 = sixT l (VR) and B AB (9) = 
exp(9(a A a B — a A a B )). The integers n A and n B repre- 
sent the number of photons subtracted, i.e. (n A , n B )=(l, 
0) or (0, 1) for the single photon detection at either site 
and (n A , n B )=(l, 1) for the coincidence detection of both 
sites. The state given by (|A2|) is a half-split state of the 
photon subtracted squeezed vacuum. Hence the situation 
described by Fig. 2£i can be equivalently transformed to 
the one in Fig. [4Jd. Although the description of (|A1[) for 
photon subtraction is only valid in the limit where the 
reflectance of the local tapping BS goes to zero, as will be 
shown in the following, the equivalency between Fig. 0^ 
and b still holds even when the photon subtraction is not 
ideal, i.e. the reflectance of the BS is not infinitesimal. 

Consider the situation in Fig. 3^, again, with defini- 
tions of modes A, B, C, and D. Here we assume that 
the reflectance of the tapping BS is finite. The following 
operator identity holds: 

B AC {9)B BD (6)B AB (Tr/A) 

= B AB (Tr/A)B AB (K/4)B AC (6)B BD {6)B AB (<K/A) (A3) 
= B AB (n/4)B AB (7r/4) 

x exp [9(a A a c - a A aJ c + a B a D - obo^)) B AB (ir/4) 



6 i ia\ ( a I - °c _ ®D , -t - ac 
B AB (ir/4) exp -a A ^ u + a\ 



- at + at „ t a c + a D 
-Ob ^ + a'r 



V2 



V2 V2 )) 

B AB (ir/4)B CD {ir/4)B AC {6)B BD {0)B% D (tt/4), 



(A4) 

(A5) 
(A6) 



where Bij stands for the beam splitter operator for modes 
i and j. Then for an arbitrary input state for mode 
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FIG. 4: a Schematic of the local photon subtraction. BS: beam splitter, HBS: half beam splitter. The letters A, B, C and D 
label the optical modes. The thinner lines represent the vacuum inputs. The gray boxes are photon detectors. Here the state 
of mode A is half-split and then photons are subtracted locally from the outputs of the HBS. b An equivalent model for a. The 
state of mode A is firstly photon-subtracted and then half-split with mode B. 



A and the vacuum states for modes B, C, and D: 

B BD (6)B AC (e)B AB (ir/4) \*) A \0) B |0) c |0) D 

= B AB {%IA)B C D (*/4)B AC (6)B bd {0)B\, d (tt/4) 
x I*)aI0) s |0) c |0) d (A7) 

= B AB (7r/4)B CD (7r/4)B AC (9) \V) A \0) B |0) c \0) D 

(A8) 

This establishes the equivalency between the models 
shown in Fig. and b again. Therefore the final out- 
put state of the protocol is identical to a half-split of the 
photon subtracted squeezed vacuum. 

Let us return to the form of (|A2|) for simplicity. Then 
we immediately have 



fit (tt/4) 1^; 



AB 



~n A +n B 

U A 



Sa(t) \0) a |0) 



B ■ 



(A9) 



which means that if we let the conditional two-mode state 
be combined at a half-beam splitter it becomes disen- 
tangled and furthermore the outputs get separated as 
a photon subtracted squeezed vacuum and the vacuum. 



Introducing new variables 



XA ± X B „ pA± PB 

x ± = ^ P± z 



V2 



V2 ' 



(A10) 



where xa,b and pa,b are the quadrature observables for 
Alice and Bob's subsystem, then from the identity (|A9[) . 
in terms of these variables the Wigner function of the 
output state can be written as 

Wout{xA,PA,XB,PB) = W v (x + ,p+)Ws(X-,P-), (All) 

where W v and W s are the Wigner functions of the vac- 
uum state and the photon subtracted squeezed vacuum 
state respectively. 

2. Local unitary equivalency between a half-split 
squeezed vacuum and a two-mode squeezed vacuum 

A two-mode Gaussian state with zero local displace- 
ments is completely specified by its covariance matrix 
given by 



(x 2 A ) \{xaPa+Paxa) \{x a x b + x B x A ) 

\{xaPa+Paxa) (Pa) \{paXb+x b Pa) 

5 (x A x B + x B x A ) \ (pax b + x B p A ) (x 2 B ) 

\z(xaPb + Pbxa) ^aPb + PbPa) \{paPb + PbPa) 



[xaPb +Pbx a )\ 



(paPb + PbPa) 
\{xbPb +Pbx b ) 
(Pi) 



(A12) 



The covariance matrix of the half-split squeezed vacuum 
|*o) = iWVWr) |0) J0) B is 



/c r cosh(r) 
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■~ r sinh(r) 



V 



o 

c r cosh(r) 
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Performing local squeezing operations on both modes, it 
can be made to have symmetric variances: 
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This is identical to the covariance matrix of a two-mode 
squeezed state with squeezing parameter |. So in terms 
of entanglement, the half-split squeezed vacuum with 
squeezing parameter r is equivalent to the two-mode 
squeezed state with squeezing parameter |. By using 
a state vector, this is described as 



S A (-r/2)S B (-r/2) |*o) = (1 - A' 2 )^ E A'™ \n) A \n) B , 

(A13) 

with A' = tanhr/2. 

3. Entanglement of the distilled states 

Entropy of entanglement is defined as the von Neu- 
mann entropy of a reduced subsystem. Every bipartite 
pure state can be brought into the form of the Schmidt 
decomposition with Schmidt coefficients c„: 



(A14) 



n=0 



From this, the entropy of the subsystem is calculated as 

E = -tv{p A \og{p A )) = -tr(p B \og{p B ) (A15) 
= -^c„logc„. (A16) 



From (|A13j) and the fact that local operations do not 
alter the amount of entanglement, we get the Schmidt 
coefficients of the half-split squeezed vacuum j^o) as 



i°) = (1 - A' 2 )A' : 



(A17) 
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FIG. 5: The entropy of entanglement for the single-photon 
subtracted states (red), the two-photon subtracted states 
(blue) and the undistilled states (black) as functions of the 
squeezing parameters r. 



Let us denote the single-photon subtracted half-split 
squeezed vacuum as Y&i). 



l* 1 > = -^a A |* >, 



sinh 



From the argument in the last section, we have 
S A (-r/2)S B (-r/2) |* x ) 

= -^S A (-r/2)S B (-r/2)a A |*„> 

= -7=7= (oa cosh(r/2) + a\ sinh(r/2)) 
x S A (-r/2)S B {-r/2) |* ) 



(A18) 
(A19) 

(A20) 



(A21) 



-^{aA cosh(r/2) + a\ sinh(r/2))(l - A' 2 )^ 
v All 



oo 



E A ' n i n >Ai»>i 



X 



(A22) 



(1 ^ E A'"(cosh(r/2)V^ |n - 1) 



+ sinh(r/2)Vn + 1 \n + 1) A ) \n) B (A23) 
(1-A' 2 )5 ' v 



nrr~ E E A ™" ' m >^ 
v-^i m=0 „=0 



(A24) 



where 



A m „ = A'"(cosh(r/2)v / n5 m ,n- 



+ sinh(r/2)Vn+ lc5 m , n+ i). (A25) 

From the singular values of matrix A mn (let them be a n ), 
the Schmidt coefficients for are given by 



(A26) 



Similar to the single photon case, we obtain a represen- 
tation of the two-photon subtracted half-split squeezed 
vacuum m the following form. 

S A (-r/2)S B (-r/2) |* 2 ) 

= -^ r S A {-r/2)S B {-r/2)a A a B |tt„) (A27) 

V 

(1 — A' 2 ) °° °° 
- ' EE^MaMb- (A28) 



iV 2 



m— n— 



where 

A/2 - 



2 sinh 4 r + cosh 2 r sinh 2 r 



(A29) 



B mn = cosh>/2)(m + l)\> m+L 5 m , n 
+ sinh 2 (r/2)mA' m - 1 c5 m ,„ 



+ cosh(r/2) sinh(r/2) x /(m + l)(m + 2)A ,m+1 <5 m+2 



+ cosh(r/2) sinh(r/2) ^(m - l)A' m_1 tf m _2, n 



(A30) 
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Then the Schmidt coefficients c„ for are given by 
the singular values (3 n of matrix B mn : 

^ = (A31) 

(ET7|) . (jAl6j) , (fMT|l and (|AT6j) we can readily calcu- 
late the entropy of entanglement for those states. Fig. [5] 
shows the entropy of entanglement for the single-photon 
subtracted states |^i), the two-photon subtracted states 
1^2) and the undistilled states | \£/ 0) as functions of 
squeezing parameter r. Note that in an actual exper- 
iment, we have several experimental imperfections and 
we end up with a mixed state output. In such cases, 
the pure state descriptions above no longer hold and the 
entropy of entanglement is not a good measure of entan- 
glement, but Fig. [5] still outlines the general behavior of 
this protocol. 

APPENDIX B: CALCULATING THE 
LOGARITHMIC NEGATIVITY 

In principle, the logarithmic negativity Ej^ can be di- 
rectly calculated when we know the 2-mode entangled 
state, as the sum of the negative eigenvalues of the par- 
tially transposed density matrix, (pab) Ta 24]. In the ex- 
periment, however, we found that Ejj is very sensitive to 
statistical noise in the measurements. Smaller data sets 
lead to larger errors in the reconstructed density matrix 
elements which ultimately leads to larger calculated £jv 
values. The intuitive understanding of this effect is the 
following: The absolute errors on each density matrix el- 
ement due to statistical measurement noise are roughly 
the same. Hence, the relative errors are large for the 
high-photon number elements which are all close to zero 
- most likely their absolute values will increase due to the 
errors. But high photon numbers contribute a significant 
amount to the overall entanglement of the state, so in the 
end more noise will give seemingly higher entanglement. 

We found empirically that the calculated negativity 
scales with the total data sample size N as Efj{N) = 
a + b/y/N. We interpret this as E^f(oo) — a being the 
"true" value that we would obtain in the asymptotic limit 



of very large data sample size, while the second term is 
the contribution from statistical noise. To obtain this 
value Ej^fipo) from a given data set of Nf u \\ samples, we 
perform multiple state reconstructions based on trunca- 
tions of the full data set. Specifically, we partition the 
full set into d subsets of Nd ~ Nf u \\/d samples each (with 
equal representation of all phase angles). The logarith- 
mic negativity is then calculated from the reconstruction 
of all d subsets, and we take the mean value of these to 
be an estimate for Ej^{Nd). We repeat the process for 
other numbers of partitions, d, and thereby get a plot as 
in Fig. [5] of the dependency of calculated entanglement 
on data sample size. A least-squares fit to a + b/\/~N 
then gives us the asymptotic estimate Em (00). We have 
confirmed by simulated data that this approach does in 
fact give the correct value for E^ within roughly ±2%. 
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FIG. 6: Estimation of the logarithmic negativity for in- 
finitely many data samples. Each point shows the average 
and standard deviation of the calculated logarithmic negativ- 
ities, Ej^{Nd) (average taken over the d subsets of the full 
data set), versus 1/y/Nd, where Nd is the size of each subset. 
The points are well fitted by a line whose y-axis intersection 
(infinite data size) gives a good estimate for the true Ej^ value 
(as confirmed by simulated data). The three plot series are 
from single-photon subtracted data sets with different initial 
squeezing levels, 2.3 (red), 3.2 (blue), and 4.2 (green) dB re- 
spectively. 
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